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Abstract 

We consider the Anomalous Hall (AH) state induced by interactions in a three-orbital per unit- 
cell model. To be specific we consider a model appropriate for the Copper-Oxide lattice to highlight 
the necessary conditions for time-reversal breaking states which are AH states and which are not. 
We compare the singularities of the wave-functions of the three-orbital model, which are related to 
the nonzero Berry curvature, and their variation with a change of gauge to those in the two-orbital 
model introduced in a seminal paper by Haldane. Explicit derivation using wave-functions rather 
than the more powerful abstract methods may provide additional physical understanding of the 
phenomena. 

PACS numbers: 
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I. INTRODUCTION 



In asking the question whether a quantized Hall effect may exist, in principle, without 
an applied magnetic field, Haldane [l| introduced an effective one-electron model on the 
two-orbital per unit-cell hexagonal crystal with complex next nearest transfer integrals. 
This has turned out to be a fecund contribution. It further augmented the topological 
arguments of Thouless and collaborators ^ for the quantum Hall effect. It also introduced 
the general discussion of the topological features of band-structures in two-dimensions and in 
particular led through further imaginative work to the suggestion and discovery of topological 
insulators {s-gI. Haldane's model also showed the connection of topological properties to 
the time-reversal- violating states due to orbital current loops in the lattice without changing 
translational invariance, see Fig. ([1]). 

Orbital current loops without changing translational symmetry were predicted to arise 
as broken symmetry states due to interactions in a three-orbital per unit cell model for 
underdoped cuprates j?-^ and have been discovered in several families of cuprates 
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12j. These loop-current states, however, do not lead to the quantized anomalous Hall effect 



(QAHE) or "Chern insulator" discussed by Haldane. The difference is that the loop-current 
states violate both time-reversal TZ and inversion X but preserve their product X. It was 



pointed out by Fradkin and Sun 



13| that there is no QAHE effect when TZI symmetry is 



present; they also pointed out that a QAHE state for the three-orbital cuprate model is 
possible, in principle through a state with symmetry different from that observed. 

In this paper, we elaborate on the work of Fradkin and Sun by deriving the Berry phase 
responsible for the QAHE state and also how the singularities of the wavefunctions vary with 
the choice of gauge. When the topologically active band (i.e., the one with nonzero Chern 
number or TKNN integer ^) is partially rather than completely filled, the Berry phase will 
still contribute to a non-quantized intrinsic anomalous Hall effect [llQ. For comparison, 
we also rederive wavefunctions for the two-orbital case. The physics of the three-orbital 
model is not surprising in the context of previous general work on topological effects in 
non-interacting lattice models; we hope that we have discussed it in a new fashion which 
will be useful to non-experts and that the particular model explains how topological effects 
might appear in three-band materials. Whenever possible we seek to explain the value of a 
Berry-phase calculation geometrically rather than simply stating the result. We also work 



2 




Figure 1: Schematic figure of loop current in hexagon lattice 

with explicit wavefunctions rather than more powerful abstract methods, although a few 
comments are provided regarding the latter. 



II. ANOMALOUS HALL STATES IN THE CU-O MODEL 

Consider the two-dimensional lattice with the structure of the copper oxides, Fig. (jj]). 
There are three orbitals per unit-cell, the ci-orbital on the copper atom and the px and Py 
orbitals on the oxygens. The minimal kinetic energy operator with a choice of gauge such 
that the d orbital is purely real and the px and the Py orbitals purely imaginary is 

Hke = it dl{sxPx,k + SyPy,k.) - t'sxSypl^^Py^k + h.c. (1) 

with Sx = sm{kx/2) and Sy = sm{ky/2) for a lattice constant taken to be 1. For simplicity, 
let the fermions be spinless. Consider only the interaction between the p orbitals. 

Hint = y^p/L^pj- (2) 

Following the procedure with which some time-reversal violating states were derived for the 
cuprates, we use the operator identity (for spin-less fermions), 

niUj = -^(Ijjjl^ + ni + nj), (3) 
3 
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Figure 2: Schematic figure of loop current in Cu-0 lattice 
where jij is the self-adjoint operator 

jij = i^clcj — CjCi). 

Decomposing^ the interaction term in ([2]), by defining 

{V/2){j,,)=zr, 

and doing a mean field calculation, one finds an additional kinetic energy term 

H'int = irc^CyPlvPy,k + h.c. 



(4) 



(5) 



(6) 



If r 7^ is a stable state, it describes loop currents flowing clockwise (or anti-clockwise) 
around the oxygen's in each unit-cell as shown in Fig ([2]). This is one of the five possible 
loop-curent states with non-overlapping loops in the Cu-0 lattice all of which preserve 
translational symmetry jl8|). In (I2]), the flux has one sign in the square formed by the 
nearest neighbor oxygens which surround a cu and another sign in the square formed by 



^ In general the diagonal in spin indices part of riirij gives (|3]), the non-diagonal can be written in terms 
of products of spin-currents with which mean-field theory gives the possibility of symmetry breaking 
topological spin-current states. The decompositions of the operator may also be done without the "i" 
in ([3]) from which various Pomeranchuk instabilities may be derived for symmetry breaking in various 
irreducible representations for spin and charge densities. 
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the nearest neighbor oxygens which do not surround a cu. As pointed out by Fradkin and 
Sun [13!] such a time-reversal violating state, which does not change translational symmetry 
or break inversion symmetry, satisfies all the conditions of a Haldane state for the Cu-0 
lattice. The other four loop-current states do not. One of those is just the photon on a 
lattice and cannot order. The other three can order and indeed order consistent with the 
symmetry of two of them (in different domains) is observed in under-doped cuprates. So, 
our consideration of states such as in ([2]) is only a specific example to illustrate the nature 
of AH states in 3-orbital models. 

We will consider the Haldane state (quantized anomalous Hall effect) of the cu-o model 
and therefore the singularities of the model with the Hamiltonian H = Hke + H'-^^. Before 
we do that, let us consider the simpler case of two-orbitals per unit-cell. 



III. TWO-BAND MODELS 



A general Hamiltonian in the space of two orbitals per unit-cell may be written as, 
ignoring an overall shift of the energy that does not affect the Berry phase and assuming 
that there is no basis of the Bravais lattice, 

[ Rs Ri-iR2\ 
H = R{k).a=[ (7) 

\Ri+ 1R2 -R3 J 
Here Ri for i = 1,2,3 are some smooth functions of and ky with period 27r. For now, 
we do not need the detailed form of these functions. It is easy to diagonalize the above 
Hamiltonian to find that there are two bands, E = ±R, with R = ^Rj + Rj + Rl 

Consider the lower band, E = — i?(k). The eigenstate can be written in two ways 
corresponding to two different choices of gauge; (the point of studying apparent consequences 
of the choice of gauge will be clear in a later section of the paper): 

They are connected by a U{1) gauge transformation 

1^^) = e"*'\^^), with e"^ = ~ f ' ^ , = - arctan(i?2/i?i) (10) 
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Then the Berry phase, is also gauge dependent, given for the two choices respectively by 



1 dRi „ dR 



At - = -;n^n^^U?2^ - (h) 



2R{R-R3)\ ^dk^ 'dk 



Af . -.(,-|V„|,-) ^ - R.^J. (12) 

y4^'s are also connected by a gauge transformation: 

Aj-< + V,,. V,^^-^(fl.||-ft§^) (13) 

The Berry curvature is gauge invariant and given by 

_dA^_dA^_J_ dRbdRc_l f^dhdRc 

~ dk, dky ~ 2i?3'-^-^«afc, dky ~ 2''""'-^'' dk, dky ^ ^ 

Here R = R/i? is a unit vector. If we integrate over the entire Brillouin zone, we find 

1 /" 72> p _ 1 f ,2, A ORbdRc 



, d kFxy = — d keahcRn-^-, — ttt- (15) 
27r J ^ 4:71 J dk^ dky ^ ' 

This is the well known result [19] that the Chern number of a two band model is 
equivalent to the winding number of the mapping from a 2D Brillouin zone which is 2D 
torus(T^) to the 2D unit sphere(S'^). This mapping can be understood by taking spher- 
ical coordinates R = i?(sin^cos0, sin^sin0, cos6'). The wave function can be written as 
10) = (— sin(^/2), cos(^/2)e*'^) which is a two component spinor. The unit sphere is just the 
Bloch sphere R = — {(j)\cr\(f)) associated with this spinor. 

The Berry phase A is a vector field defined on the momentum space or the torus. Since 
the torus is mapped to a sphere, one can also think of the Berry phase as defined on this 
sphere. In spherical coordinates, we have 

A^ = 1±^0, A^ = -1^^0 (16) 
2sin^ ^' 2sin^ ^ ^ ' 



This is just the vector potential of the Wu-Yang monopole 20| on a unit sphere. The 
magnetic field is B = V x A = ■^^^{A^sm.6)v = — |r. The magnetic field is like that of 
a monopole with charge g = —2tt. Furthermore, this implies that a monopole-like singular 
point is located at the center of the sphere Ri = R2 = R3 = 0, which is also the point at 
which the two bands become degenerate. 

The winding number is easier to compute than the Chern number, since it is directly 
written in terms of matrix elements of the Hamiltonian and does not require computing the 



eigenvectors. If Ri, R2 and -R3 are independent of each other, and each of them can take 
both positive and negative values at the point where other two components are zero, then 
R will sweep out the whole unit sphere which encloses the singular point inside it. Then the 
winding number or the Chern number is nontrivial. But in general it is also possible that 
R winds the sphere twice in opposite directions and cancel out or comes back to a given 
point as k is varied over all its values without sweeping the entire sphere. Then the Chern 
number is zero. 



As a specific example 2l|], take R = (sin k^, sin ky, m + cos kx + cos ky). For < m < 2, 
there is always a gap between the two bands and it is easy to see that R sweeps the whole 
sphere. One can also directly verify that 



1 [ i2i "^os + COS ky + m cos kx cos k 



y 



-1 



(17) 

[sin^ kx + sin^ ky + {m + cos kx + cos kyy]^/'^ 

Consider next the Haldane model It has a staggered flux inside each unit cell. The 
Hamiltonian is 

H = ti(^Y^ cos(k ■ ai)a^ + ^ sin(k ■ ai)o-^ j + (^M - 2t2 sin (p ^ sin(k ■ hi)^a^ (18) 

i i i 

The maximum of ^■sin(k ■ hi) is 3-\/3/2, thus for \M/t2\ < 3-\/3| sin0|, R will cover the 
whole sphere. Indeed, for > 0, one can directly verify that c = 1 in this case. 

The condition on M is simply the necessary condition for the AH state that the monopole 
singularity exist which requires that -R3(k) go through zero at some point k and change sign 
as that point is crossed. 

Fradkin and Sun [13] pointed out that if both time-reversal TZ and inversion X are broken 
but the product TZX conserved, there can be no anomalous Hall state: the diagonal compo- 
nents of the Berry vector potential must vanish. Note that TZX Rsa^ = —R^cr^. Therefore 
if TZX is conserved -R3(k) = for all k. Then the sphere -R(k) turns to a circle and no 
singularity can be defined. 

In contrast to this is the 611 type loop current states realized in the three-orbital model 



for cuprates 



IOI-I12II. Its mean field Hamiltonian in the same basis as (H]) is 
/ 

H = 







ttSx I Cx '^t Sy I l/T" Cy 

—itSx — irCx 



\ 



v 



S y liir Cy 



t SxSy 



t SxSy 





(19) 



If we define id^ — d\ then tfie Hamiltonian in the new basis is a real matrix. Actually, if 
TZI is invariant, then the Hamiltonian (in momentum space) can always be written as a 
real matrix and the phase of all eigenvectors are constant. Therefore such a state though 
violating time-reversal cannot have an AH state. 



IV. THREE BAND LOOP CURRENT MODEL WITH AH 



A. Chern number and winding number 

Now we turn back to the three-orbital copper-oxygen model given by Hamiltonian H 
Hre + H' , which in the space of d,Px,Py is 



H 




-itSj 







- ITCxCy 



t SxSy ~\~ iVCxCy 





(20) 



with Sx — sm{kx/2) and Cx — cos{kx/2), etc. Note that even though the matrix elements as 
a function of kx and ky do not have the period 27r, the energy dispersion as a function of k^ 
and ky does have a period 27r. As commented on below, the Hamiltonian for a tight- binding 
model with a basis is not strictly periodic when k is translated by a reciprocal lattice vector 
G, but rather is transformed by the unitary matrix exp(iG • a^, where is the location of 
the ith site in the unit cell. 

To diagonalize the above Hamiltonian, one has to solve a cubic equation to find out 
the eigenvalues. Since Chern number is topological invariant, if we deform the Hamilto- 
nian without the bands crossing, the Chern number will stay the same. We can therefore 
simplify the problem by dropping the t'sxSy term and come back later to ensure that this 
simplification is valid. The simplified Hamiltonian is 



H' 



( 



\ 





—its,. 



itS-r 



\ 



'tVCx Cy 







(21) 



This Hamiltonian can written as H — R ■ h with Ri — tSx, R2 — —tSy, Rs — rCxCy and 
Lx, Ly, Lz are the spin 1 representation of the SU{2) generators in contrast to the spin 1/2 
representation of the SU (2) generators for the two band case. 
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We will discuss the location of singular points in simple gauges for model H' of Eq (4). 
They should be qualitatively similar to model H of Eq (3). It is easy to find the eigenvalues, 



£J = 0, ±i? with R = ^JR\ + Rl-\- Let us focus on the lowest band Ei = -i?, the 
corresponding eigenvector can be written in two different ways labeled by A, 5, 

1^5") = ^2^o - ^3, -R^R^ + ^^^1, -^1^3 - iRR2Y (22) 

^J2[R^ - Rl)R ^ / 

IVf ) = , ( - RiR^ + iRR2, -R1R2 - iRRs, R^ - RlY (23) 

y 2(i?^ — R()R ^ ^ 

These two wave functions are related by a C/(l) gauge transformation 

iV^n^e^l^^), ,^<t> = ^M^±^£^=, = -arctanf|^) (24) 
' ^ ' ^{R^-Rl){R^-Rl) \R^rJ ^ ' 

The Berry phase is given by 

They are also connected by a U{1) gauge transformation, 

+ (27) 

The Berry curvature is given by 

Note there is a factor of 2 difference compared to the two-band model. The Chern number 
is the integral of the Berry curvature in the BZ 

\ f \ f OR OR 

27ry 4[(is,)2 + (is,)2 + (rc,c,)2]3/2 ^ ^ 

In this model, -R3 = rc^Cy takes only positive values and Rx = ts^ and R2 = —tSy continu- 
ously vary from —t to t and do not go back. Thus R only sweeps half of the unit sphere. 
Thus the surface integral J (fk'K{QxK x 5j,R) = 27r and Chern number c = —1. 
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We have obtained the relation between the Chern number and winding number for the 
spin 1/2 and 1 representations of SU{2) algebra. Actually, we can generalize this relation 
to any spin n/2 representation of SU{2) for integer n. The Hamiltonian is given by 

H = RiJx + R2-Jy + Rs-Jz (30) 

Here Ja for a = x,y,z are nhj n matrices and satisfy [J^, Jb] = i^abcJc- Then the eigenvalues 
are Ei/R = —n/2, —{n/2 — 1), ■ ■ ■ ,n/2 — l,n/2. Then for the ith band, we have 



1 



2tx J ~' '"'""^ R 27r 
To prove the above result, it is better to consider H/ R = R ■ J which has the same Chern 
number as H. Since R is a unit vector, it can be parameterized by spherical coordinates 



d kdahcRa 



dRb OR, 



(31) 



R = (sin cos 0, sin ^ sin 0, cos 6'). Here 6, are functions of and ky. 

To be specific, we take Jz as a diagonal matrix = diag{n/2, (n/2 — 1), 
1, —n/2}. We can make a rotation to diagonalize H/R as. 

For the ith band, we have the eigenvector 

= e-''t''^e-'''yYi,, ni(0, • ■ • , 1, • 



-n/2 - 



(32) 



(33) 



In the above vector, only ith component is 1 and all others are zero. The Chern number 
can also be written in terms of differential forms as 



c = / d k. ^ „ 

27r / V okr ok 



dip'^ dip dip^ dip \ 



dky dkx / 27r 



dip^ A dip 



(34) 



It is easy to find 

dip = (^-iJze-'^^^e-'^^yd(P-ie-''^-^^e-'^^yJyde^ni 
dip^ = nj{ie'^^ye'^^'Jzd(p + iJye'^-^ye'^-^^de^ 
Using the above in Eq (IMll . we find 



c = -— / (jyc'^-'y Jzc-'^^ydd Ad(P + e'^^y JzC-'^-'y J, 



/ Tr 

27r 



y 

idj; 



A de n,- 



6 ^^JyJz JzJy)^ 



-''^ ■ (n.nf ) 



dO A d(p 



1 

2^ 
1 

2^ 



Tr 
Tr 



(n^nf) 



[cosOJx + sin^Jj 



d9A 



dO A d(p 



(35) 
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The matrix (rijnf ) has 1 as its ith diagonal element, and all other elements are zero. Thus 
for any nhj n matrix A, we have Tr[y4(njnJ')] = An. Since only has off-diagonal elements, 
only the second term of Eq (l35i) contributes. Therefore, 



1 /" E- 1 C 

c = — {J^)ii / sin 6'd6' Ad(j)= — — / sin 6'rf6' A d(f) (36) 

On the other hand, the winding number integral in terms of spherical coordinates can be 
written as 

/dR OR 1 /" /" 
d'^keabcRa-Q^-Q^ = 2 J ^^b'^^^^^b ^ ^-^c = j ^^"^^dO A d(j) (37) 

Combining the above two equations, we find the desired results. 

There is a subtle point about this loop current model's behavior when k is shifted by a 
lattice vector. Recall that Ri = tsm{kx/2), R2 = — t sin(/cy/2), R^ = r cos{kx/2) cos{ky/2) 
and they are functions of kx and ky with period An not 27r. Thus the Hamiltonian is not 
invariant under changes k^ ^ k^ + 2v: and ky — ky ^ 2Ti . This happens generally if the 
lattice contains a basis. It is easy to see that the eigenstates are also not invariant under 27r 
shift. But if we rewrite the wave function in real space, we find 



(38) 



-(I sin A;^ - itR) sin(A;^./2)e^'=-/2 
\^ -(f sin kx + itR) sin(A;j,/2)e^^^/2 y 

Since R and -^1^2,3 has period 27r, the above wave function is invariant under 27r shift as one 
expected. 

If we define the torus to be — vr < k^^y < vr, then R^ is always positive and the torus 
is mapped to half of a unit sphere by R. It might seem that in this case, the monopole 
singular point is not enclosed by the half surface, therefore the Chern number need not to 
be quantized and may have no topological meaning. But if we look at the vector potential 
in Eq.( l26|) . they are invariant under 27r shift, and thus are well defined on the whole torus. 
Therefore, according to Dirac's arguments, the integral of field strength on this closed surface 
should give a quantized topological invariant. In the mapping Il{kx,ky), the boundary of 
BZ is mapped to the equator of the sphere. As discussed above, the wave functions are the 
same at the two points like (— vr, ky) and (vr, ky) on the boundary of BZ. These two points 
are mapped to {—T^,Sy,0) and {tt-^, Sy,0) on the equator. These two points should be 
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identified, and then the boundary of the half sphere is glued together to make a closed 
surface which is topologically equivalent to a sphere. In this loop current model, we have 
wound over this closed surface once. 

B. Singular points of wavefunctions and their gauge-dependence 

We now discuss the singular points of wave function inside the BZ for the three-orbital 
case. One can see that -0^^ is well defined for all possible values of -Ri,2,3 except when 
i?3 = ±R or Ri — R2 — 0. So this is the singular point of ipf. It corresponds to — ky — 
in the BZ. Similarly, t/jf has singular points when Ri — ±R or R2 — R^ — 0. It correspond 
to kx = ±TT, ky = in the BZ. Therefore the location of singular points depends on the 
choice of gauge. 

The reason that any gauge has singular points is simply that the Chern number being 
nonzero implies no continuous gauge can cover the whole Brillouin zone. It is easiest to 
understand the connection of these two by thinking of a spherical rather than toroidal 
Brillouin zone. If a single gauge covered the whole sphere, then we could apply Stokes' 
theorem to relate the Chern number, which is the integral of the Berry curvature over the 
whole sphere, to the integral of the Berry connection around a tiny circle, which must be 
zero. In the same way, a nonzero Chern number is an "obstruction" to having continuously 
defined wavefunctions over the whole Brillouin zone. 

Each form of wave function is valid only on one patch of the torus and the two are 
connected by C/(l) gauge transformation on the boundaries. This is exactly the same as for 
the Wu-Yang monopole. ip^ and ipf define a U{1) bundle on the torus. Let be the open 
set which covers the torus without the point k^ = ky = 0. Let be the open set which 
covers the torus without the point k^ = ±71, ky = 0. Then ip^ is well defined on U"^ and 
■0f is well defined on U^. Both and ip^ are well defined on the overlap U^f]U^ and 
they are related by C/(l) gauge transformation — e^'^ipi . We can take a closed loop as 
the boundary of C/^ and C/^, then we find a mapping from this closed loop to U{1). Then 
the Chern number is just the winding number of this mapping. 

We can make this connection more explicit by taking a small closed loop as c k'^ -\- ky — 6^ 
with small positive e to enclose the singular point of ip^. Making a small k expansion, we 
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find i?i,2 ~ ^x,?;/2, R3 ^ 1, R ^ 1. Then we have 

'^--^r^tan^, V,0^(^^,-^A^) (39) 
Thus the Chern number can also be obtained as a line integral around this small loop 

Even without detailed calculation, on can see that as the momentum goes around this loop, 
the phase difference (p of Eq flMj) also goes around one circle. Thus the winding number is 
— 1, which agrees with the previous calculation. 

The Chern number can also be defined as a line integral c = ^ ^^Afj_dl^, where integral 
path C is the boundary of the BZ. Since the Berry vector potential depends on the gauge 
choice, to get a correct answer for the Chern number, one has to choose a gauge such that 
the wave function has no singular point inside the loop. In the Cu-0 model, we should use 
A^. It is easy to see that = for ky = ivr and — tt < < it and Ay = for k^ = ivr 
and — TT < ky < 7T. Thus we find loop integral f Andlfj_ = 0. But since there are two singular 
points kx = ±vr, fc^, = on the boundary, we should choose an integral path with a small 
semi circle to circumvent the two singular points. Expanding Aj^ around kx = ±7r, ky = 0, 
we find 

Aj(±7r + kx, ky) = (^kx, -ky^ (41) 

Let e denote the small circle around the singular point, then 

^ I A,dl, = I kxdk - k dkx ^ _^ ^^2) 



27r/c " 27^1 kl + kl 

This result agrees with the Chern number calculated from the Berry curvature. 

For the middle band E = and the wave function is |?/'2) = ^{^3^ -^2? -Ri) which is always 
real in the BZ. Clearly this band is topologically trivial and Chern number c = 0. The sum 
of the Chern numbers of all bands is always zero. Thus for the top band we have c = 1. 

Now we come back to the original Hamiltonian H. We show that if t' is not too large, 
there is no band crossing so the result stays the same as H'. Let i?4 = ts^Sy, the eigenvalues 
are the roots of cubic equation 

- {Rl + R^)E + 2R4R1R2 = (43) 
13 



The condition for degeneracy is — [(-R4 + -R^)/3]^ + {RiRiR2Y = 0, It can be rewritten as 
-[{Rl + l)/3f + (i?4i?i4)' < -[{Rl + + Rl/A = 

with R4 = R4^/R. The last equahty of the above is satisfied only for Rl = 1/2. In our model, 
for simplicity, we take r = t, then the maximum value of R4 is ^ ai = ky = tc. Thus as 
long as t < t', there is no band crossing and the Chern numbers stay the same as H'. 

It follows from the discussion above that if the band is only partially filled, i.e. there 
is a Fermi surface, there can be no topologically protected currents. Nor can there be 
any singular properties of the Fermi liquid coming from the physics of the (non-quantized) 
anomalous Hall effect because one can always move the singularities of the wave-functions 
away from the Fermi surface by a suitable choice of gauge. The intrinsic contribution to the 
anomalous Hall effect can be calculated by integrating the Berry curvature over the Fermi 
surface volume. This procedure is best carried out numerically by writing the Chern number 



in terms of the projection operator [22|] onto the occupied subspace, Pk = \4'k){4'k\, as this 
object is manifestly gauge-invariant. 

V. CONCLUSIONS AND REMARKS ON THREE-DIMENSIONAL MATERIALS 

We have discussed the calculation of Berry phases and the contribution to the Hall effect 
in a three-band model motivated by the copper-oxygen planes of the high-temperature su- 
perconductors. There are two topologically active bands (i.e., with nonzero Chern number 
±1) separated by a middle band with Chern number zero. We have neglected interaction 
effects, but there is considerable recent interest in the possibility of fractional quantum Hall 
phases when a band of nonzero Chern number is partly occupied. Bands with Chern number 
±2 or larger are particularly interesting as, while these are formally equivalent to bilayer 
quantum Hall systems, it has been argued that they are likely to support novel fractional 



quantum Hall states because the nature of interactions is modified [23[. That work proposed 
creating the Chern number ±2 bands via an oxide heterostructure, and this technique might 
also enable the PT-breaking model described here to be realized at the interface between a 
cuprate and another material. 

In closing, it may be useful to mention another topological property enabled when a model 
has three bands rather than two; the CP violation enabled by the three-by-three Cabibbo- 
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Kobayashi-Maskawa matrix in the Standard Model provides a well-known example of how 
three-by-three matrices allow additional subtleties. The orbital contribution to the linear 
magnetoelectric effect, i.e., the polarization induced by a magnetic field dPi/dBj or magne- 



tization by an electric field dMi/dEj, has recently been a subject of active study 2l|, |2J-|26|. 
The topological part of this effect is the scalar diagonal part ("axion electrodynamics"), 
computed by the Chern-Simons integral over the Berry connection 2l| , which vanishes in a 
purely 2D model such as that considered here. 

In a 3D model, either P or T symmetry quantizes the scalar diagonal part to only 
two possible values, corresponding to ordinary and topological insulators. Without these 
symmetries, in a 3D model with only two bands, the Chern-Simons integral is quantized and 
computes the "Hopf invariant" of the band structure, viewed as a mapping to the sphere . 
Three bands are required in order to generate generic values of the scalar magnetoelectric 
coupling. The computation for a single occupied band, in which case the Chern-Simons 
integral is Abelian, has been recently discussed [28|. We hope that the results of this 
paper lead to further study of the consequences of wavefunction topology for transport and 
magnetism of oxide materials. 
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